Bond-bending modes and stability of tetrahedral semiconductors under high pressure: 

a puzzle of A1N 
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Lattice vibrations of the wiirtzite-type A1N have been studied by Raman spectroscopy under high 
pressure up to the structural phase transition at 20 GPa. We have shown that the widely debated 
bond-bending E\ mode of w-AIN has an abnormal positive pressure shift up to the threshold of the 
phase transition, whereas in many tetrahedral semiconductors the bond-bending modes soften on 
compression. This finding disagrees with the results of ab initio calculations, which give a "normal" 
negative pressure shift. Combination of high dynamical and low thermodynamical stability of A1N 
breaks the correlation between the mode Griineisen parameters for the bond-bending modes and 
the transition pressure, which holds for CdS, InP, ZnO, ZnTe, ZnSe, ZnS, Ge, Si, GaP, GaN, SiC 
and BeO. 

PACS numbers: 



For the last four decades, a rather curious and at first 
sight surprising phenomenon of softening of tetrahedral 
semiconductors under compression has been reported in 
a number of experimentalllUl EJ1J1 II 0J1 IS 
and theoretical [U [H HOl ElM Hill works. 
Thus, Si, Ge, A n iB v and A n B V i semiconductors have 
pressure-sensitive phonon modes with negative pressure 
shifts [3J, \A lj- liil • These soft modes are shearing modes, 
involving bond bending in the first order of the strain 
|l5l Il8l|. Phonon frequency drop, more pronounced for 
high-Z materials, reaches about 30% at the threshold 
of the pressure induced phase transitions, when covalent 
tetrahedral structures lose their stability and transform 
to more densely packed arrangements. 

Earlier, the only known example of a bond-bending 
mode having a positive, although weak, pressure shift 
was being the "pure bread" bond-bending E\ mode of 
the low-Z wiirtzite-type BeO 0. This result was being 
regarded as an exception until recent Raman measure- 
ments found a similar behavior for the E\ mode of the 
wiirtzite-type A1N [2(j ■ Ab initio calculations (2jJ, 0] 
did not catch this feature, giving a "normal" negative 
pressure shift for the E\ mode of A1N. This discrep- 
ancy deserves special attention, since the bond-bending 
elasticity is one of the most prominent manifestations of 
directional covalent bonding, and its pressure behavior 
should be well accounted for. The point to be made here 
is that A1N represents a complicated case of covalent ver- 
sus ionic bonding p3 | : although its valence charge distri- 
bution is highly io nic [2^, A1N adopts the tetrahedrally 
coordinated wurtzite structure under ambient conditions 
|24|. Aspiring to clarify the question, we have under- 
taken a complementary high-pressure Raman study of 
vibrational modes of the wiirtzite-type A1N up to its sta- 
bility limit at about 20 GPa. The pressure dependence 
of the low-frequency bond-bending E\ mode is traced up 
to the threshold of the pressure induced phase transition 
for the first time. 

The A1N samples were 20 /xm-thick crystals grown on 



the sapphire substrate by vapor phase epitaxy. Pres- 
sure was produced using the diamond-anvil pressure cell. 
Compressed helium and methanol-ethanol mixture were 
used as a pressure-transmitting medium in the first and 
in the second experimental run, respectively. Pressure 
was measured in situ by the ruby luminescence technique. 
The Raman spectra were measured using the THR-1000 
triple spectrometer equipped with an OSMA detector 
(the first run), and the Dilor XY double spectrometer 
equipped with the CCD detector (the second run). An 
Ar + laser (A = 514.5 nm) was used as a source of exci- 
tation. All spectra were recorded in the backscattering 
geometry at ambient temperature. 

For the hexagonal wurtzite structure with the space 
group P63mc {Z = 2), a factor-group analysis predicts 
the following six sets of optical modes at the k = |25| : 

T op = A 1 +2B 1 + E 1 + 2E 2 , 

where A\, E\ and Ei are Raman active modes, and B\ 
modes are silent. A\ and E± are also infrared active, 
and split into the longitudinal and transverse components 
(LO and TO). The lowest-frequency E\ mode is a bond- 
bending mode. 

The Raman spectrum of w-AIN has been measured 
previously under ambient conditions and analyzed in 
some detail, inclu ding the effects of polarization and 
anisotropy (26L 1271 |28|. Our ambient pressure Raman 
frequencies v E \ — 249 cm -1 , va 1 {to) — 610 cm -1 , 
v E 2 = 657 cm -1 , v El (TO) = 669 cm" 1 , va^lo) = 890 
cm -1 , and v E „<T t p\ = 910 cm -1 agree with very reli- 
able data of Ref . |2rj I27L l28l | within 1%. On increase in 
pressure, all Raman bands shift continuously to higher 
phonon energy with no intensity loss to about 18 GPa. 
Above 18 GPa, the bands weaken and disappear at about 
21 GPa in both experimental runs due to the phase tran- 
sition to the rock salt structure |29j, |30j . Representative 
Raman spectra of w-AIN in the low-energy region as a 
function of pressure are shown in Fig. ^ 
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FIG. 1: Raman spectra of AIN as a function of pressure in 
the low-energy region. The spectral resolution is 0.5 cm -1 . 



Fig. compares measured and calculated pressure de- 
pendences of the Raman frequency for the E\ mode. The 
mode pressure coefficients v\ calculated using the linear 
least-squares fit vi = z/? + v[P, where is the frequency 
of the mode i at pressure P, are listed in Table I |31j . 
The pressure dependence of the E\ frequency, measured 
in our experiment, is weak but apparently positive and 
linear at pressures to the threshold of the phase transi- 
tion [32. This finding confirms the results obtained in 
[2fj|. although our pressure slope of the E\ frequency is 
somewhat lower (see Fig. EJ. So, the positive pressure 
shift of the E\ frequency of AIN should be regarded as a 
reliable experimental fact. 

As we have mentioned above, ab initio calculations 
do not reproduce this feature of AIN, although 
the agreement between measured and calculated pressure 
shifts for the high-frequency bond-stretching TO and 
LO modes is very good (see Table I). Certainly, the value 
of the discussed effect for the E\ mode is quite small, and 
the lack of ab initio calculations is not so crucial. Nev- 
ertheless, this lack is very instructive to realize a subtle 
pressure-sensitive balance between stabilizing and desta- 
bilizing contributions to the bond-bending force constant 
I l'2t Il4l Il5| and may serve a further improvement of the 
theory. 

Table II compiles the commonly used mode-Griineisen 
parameters ji — — din Vi/d In V obtained for the bond- 
bending modes i in a scries of tetrahedral compounds 
[33|]. Negative value for 7, is observed in each case, ex- 
cept for SiC, BeO and AIN. Although SiC has an essen- 
tially zero 7^, its quadratic pressure coefficient is neg- 
ative Q. Thus, BeO and AIN appear to be the most 
stable materials with respect to the bond-bending mode 



FIG. 2: Comparison between the measured and calculated 
pressure dependence of the Raman frequency for the E\ mode. 
The squares are the first run data, obtained with compressed 
helium as a pressure transmitting medium. The triangles are 
the second run data, obtained with ethanol-methanol mixture 
as a pressure transmitting medium. The solid line 1 is a linear 
fit of the first run data. The dotted line 2 is guided for the 
eye. The lines 3 and 4 are calculated dependences obtained 
in [2l| and [2fj|. respectively. The line 5 is the experimental 
dependence obtained in [2(l |. All data are shifted along the 
vertical axis in order to coincide the 1 atm frequencies with 
the value 249 cm" 1 obtained in our first experimental run. 



TABLE I: Measured and calculated linear pressure coefficients 
v[ (cm -1 GPa -1 ) for the AIN Raman frequencies. 

E\ A x {TO) El Ex{TO) Ai(LO) E-^LO) 



Expt. 

Run 1 
Run 2 
Gofii a 
Calc. 

Gorczyca 6 
Goni° 



0.05(1) 
0.05(1) 
0.12(5) 

-0.29 
-0.03 



3.8(2) 
4.3(2) 
4.4(1) 

4.29 
3.0 



4.9(2) 
4.65(3) 
4.99(3) 

4.79 
4.2 



4.5(1) 
4.55(6) 
4.55(3) 

4.36 
3.8 



4.0(1) 



3.5 



3.6(7) 
4.61(1) 



4.0 



"Rcf. 20 
'Ref.m 



on compression. 

This bond-bending stability of AIN looks rather sur- 
prising, since its wiirtzite phase has a rather limited sta- 
bility range up to 21 GPa. Indeed, experimentally it has 
been found that the stability of tetrahedral structures 
with respect to the bond-bending modes correlates with 
their absolute stability under pressure in such a way that 
the frequency drop is faster for less stable compounds. 
Weinstein 0, Q has discovered that for six diamond and 
zinc-blende structure semiconductors ZnTe, Ge, Si, ZnSe, 
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TABLE II: Mode-Griineisen parameters for the bond-bending 
modes. 



Material 


Structure 


Bond-bending mode 


Ji 


CdS 


wiirtzite 


E 2 


-2.7" 


InP 


zinc blende 


TA(L) 


-2.0" 


ZnO 


wiirtzite 


E 2 


-1.8" 


GaAs 


zinc blende 


TA(L) 


-1.7 


ZnTe 


zinc blende 


TA(L) 


-1.5 


Ge 


diamond 


TA(L) 


-1.52 6 


Si 


diamond 


TA(L) 


-1.3 C 


ZnS 


zinc blende 


TA(L) 


-1.18" 


GaP 


zinc blende 


TA(L) 


-0.81 d 


GaN 


wiirtzite 


E 2 


-0.426 6 


SiC 


hex.(6H) 


E 2 


0.0 / 


BeO 


wiirtzite 


E 2 


0.04 9 


A1N 


wiirtzite 


E 2 


0.04' 1 



"Ref. 
<>Ref. 
c Ref. 
d Ref. 
e Ref. 
•^Ref. „ 
9 Ref.|3 
''This study. 



ZnS and GaP there is a remarkable linearity between the 
mode Griineisen parameter jta(x) for the bond-bending 
TA(X) mode and the transition pressure Pt r for these ma- 
terials. Ab initio calculations by Yin and Cohen have 
given some insight into the nature of this correlation. 
They have shown that those individual contributions to 
the total crystal energy of Si and Ge, which stabilize 
the bond-bending TA(X) mode, are just the same ones, 
which preserve the diamond structure from rearrange- 
ment to more densely packed structures 0- The bal- 
ance between stabilizing and destabilizing contributions 
drastically depends on the specific volume, destabilizing 
contributions becoming stronger under compression. 

At present it is possible, following Wcinstein and Zallen 
0, , to trace the 7$ — P tr correlation to much higher 
pressures, using recent data for w-GaN j8j, SiC-6-ff 0, 
l3l| and w-BeO [3(|. The mode Griineisen parameters 
determined for the bond-bending TA(L) phonons of InP, 
GaAs, ZnTe, ZnS, Ge, Si, GaP and for the bond-bending 
E 2 phonons of CdS, ZnO, A1N, GaN, SiC and BeO are 
plotted in Fig. [3] as a function of the transition pressure 
P tr ranging to 140 GPa. 

Fig. |3 shows that the 7, versus Pt r points form a 
smooth curve deviating from linearity at P tr > 20 GPa 
and approaching to zero at P tr w 100 GPa. On the basis 
of this curve one may reasonably expect that materials 



with ji > should be at least as stable under pressure 
as SiC does. However, the point for A1N falls far afield 
from the common curve and tends to continue a linear 
dependence drawn by high-Z materials. 

So, we see that a simple universal correlation between 
the pretransitional behavior of the bond-bending modes 
and the stability of tetrahedral semiconductors under 
0,5 
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FIG. 3: Correlation between the mode Griineisen parameter 
7i for the bond-bending TA(L) [E 2 ] modes and the transition 
pressure in a series of tetrahedral semiconductors. Solid and 
dashed lines are guided for the eye. References are given in 
Table II and in the text. 



high pressure does not exist. To our knowledge, the ques- 
tion why AIN combines a high dynamical and a low ther- 
modynamical stability remains still open, and a detailed 
microscopic treatment, possibly along the direction indi- 
cated by Yin and Cohen |15| , is required. 

The authors wish to thank A. Dobrynin for growing 
the AIN crystals. E. V. Iakovenko is grateful to A. F. 
Goncharov for his assistance in Raman measurements. 
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